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ABSTRACT 

We report on the first ground-based measurement of the relativistic beaming 
effect (aka Doppler boosting). We observed the beaming effect in the detached, 
non-interacting eclipsing double white dwarf (WD) binary NLTT 11748. Our 
observations were motivated by the system's high mass ratio and low luminos- 
ity ratio, leading to a large beaming- induced variability amplitude at the orbital 
period of 5.6 hr. We observed the system during 3 nights at the 2.0m Faulkes 
Telescope North with the SDSS-g' filter, and fitted the data simultaneously for the 
beaming, ellipsoidal and reflection effects. Our fitted relative beaming amplitude 
is (3.0 ± 0.4) x 10~ 3 , consistent with the expected amplitude from a blackbody 
spectrum given the photometric primary radial velocity amplitude and effective 
temperature. This result is a first step in testing the relation between the pho- 
tometric beaming amplitude and the spectroscopic radial velocity amplitude in 
NLTT 11748 and similar systems. We did not identify any variability due to 
the ellipsoidal or reflection effects, consistent with their expected undetectable 
amplitude for this system. Low-mass, helium-core WDs are expected to reside 
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in binary systems where in some of those systems the binary companion is a 
faint C/O WD and the two stars are detached and non-interacting, as in the case 
of NLTT 11748. The beaming effect can be used to search for the faint binary 
companion in those systems using wide-band photometry. 

Subject headings: stars: individual (NLTT 11748) — white dwarfs 



1. Introduction 



The movement of a light-emitting object relative to an observer causes relativistic beam- 
ing of the light (aka Doppler boosting) since the specific intensity is not a relativistic in- 
variant. When the object's radial velocity (RV) is periodically modulated, as is the case in 
binary systems, a periodic sinusoidal variation in the observed flux is induced. This effect 
is one of three that can be seen in light curves of binary systems, not necessarily eclips- 
ing. The other two are the ellipsoidal effect, when a member of the binary system is tidally 
affected by the gravitational pull of its companion, and the reflection effect, when light orig- 
inating from one component is re flected by the other. The three effects were discussed in 
th e context of planet-h ost stars by lLoeb fc Gaudil (120031 ). and main sequence stellar binaries 
by IZucker et al.l ( 120071 ) . Both studies addressed binary systems with orbital periods much 
longer than a day, finding that the required photometric accuracy necessitates space-based 
observations. 



Observational detections of the beaming effect are rare. iMaxted et al.l ( 120001 ) noted that 
the beaming effect is likely to be seen in their ground-based data of a sdB star + WD binary 
system with a period of P=2.28 hr, although the effe ct's amplitude was not me asured and 
the light curve is dominated by the ellipsoidal effect. Ivan Kerkwijk et al.l (120101 ) measured 
the beaming and ellipsoidal effe cts in a Kepler lig ht curve of an A-type star + low-mass WD 
binary system with P=5.19 d (IRowe et al.l 120101 ), and in fact used thes e measurements to 
resolve the puzzling nature of the system (see also lEhrenreich et al.ll2010l ). iMazeh &: Faigler 
( 120101 ) identified the beaming a nd ellipsoidal e f fects in CoRo T data of a brown dwarf host 
star (CoRoT-3; P=4.26 d), and lBloemen et al.l ( 120101 ) observed all three effects with Kepler 
in a 9.69 hr sdB + WD binary. 

We show here that for high mass ratio detached, non-interacting double white dwarf 
(WD) binaries the relative beaming amplitude is above the 10 -3 level, while the ellipsoidal 
and reflection amplitudes are smaller by an order of magnitude or more. Therefore, these 
systems represent an excellent opportunity to measure the beaming effect using ground-based 
telescopes. 
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Following; the discov ery of the first eclipsing detached double WD binary NLTT 11748 
by ISteinfadt et al.l (j2010[ ) we immediately noticed it is an excellent candidate for ground- 
base doj3serva_tioj^ is large (~270 km 



-l. 



Steinfadt et al.lboiol : ETawka et alJlioiot iKilic et al.ll2010bl ). while the secondary is 



^30 

times fainter than the primary Jsteinfadt et alJboioh . In addition, the orbital period is short 
enough to be covered within a single night. 

NLTT 11748 and similar systems where both the (photometric) beaming amplitude and 
the (spectroscopic) RV amplitude can be measured are a valuable tool for testing the relation 
betw een these two quantities, which is currently ap proximated using theoretical assumptions 
(e.g., 



Zucker et al.l 120071 : Ivan Kerkwiik et al.ll2010f ). 



We discuss the beaming effect in non-interacting double WD systems in § [2J Our ob- 
servations and data analysis are described in § [3J We present the results in § 0] where we 
also test them and compare our findings to theoretical expectations. In § [5] we summarize 
our results and give our conclusions. 



2. Beaming in non-interacting double WDs 



The relatively large beaming amplitude in non-interacting double WD binaries originates 
from the unique WD mass-radius, or mass-luminosity relation. In those binary systems the 
photometric primary (hereafter, the primary) is a helium core WD, a few times smaller in 
mass but larger in radius and luminosity than t he photometr ic secondary ( hereafter, the 



second ary) , a C/O core WD (see also Table 1 in IKaplanl 120101 ) . As shown by IZucker et al. 



(120071 ) the observed beaming effect is the weighted difference between the individual beaming 
effects of the two binary components, shifted by a phase of 0.5. However, the larger mass 
and smaller luminosity of the C/O WD make its own beaming effect much smaller, and even 
negligible relative to that of the He WD. Th is can be seen by a ssuming both stars radiate 
as blackbodies and rewriting Equation (6) of IZucker et al.l (120071 ) as: 



A 



beaming ^1 



a 2 mi Fvp 
a[ m 2 F v ,i 



c 1 + F V JF V - 



where K\ is the radial velocity amplitude of the primary, c the speed of light, m\ (m 2 ) is the 
mass of the primary (secondary) and F Vj i {F v ,2) is the flux from the primary (secondary). 
a[ 2 relates to the spectral index 0^2 = 



rflog(F„,i, 2 ) 
dlog(u) 



as: 



a 



1,2 



«1,2 
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where x^2 = hv jhT^^^ and T^i (T^^) is the effective temperature of the primary (sec- 
ondary) . 

For a high mass ratio non-interacting double WD binary the term 2 /f"\ m Eq. ([I]) 
right hand side becomes of order unity. This follows from the mass ratio and the flux ratio 
of the two stars, and also from the fact that oi 2 /a x is itself close to unity for the expected 
effective temperatures of the two stars. Therefore Eq. (jTJ) reduces to: 

A(u) beaming « a[ K± = 10- 3 a[ ^ ^ ^ , (3) 
or in physical parameters: 

A (v) ^ m , u 2 x 10-' a[ £ (=«+2S) V3 (j^y V3 -»« , (4) 

where i the orbital inclination angle. The approximation neglects the contribution of the 
secondary and assumes a blackbody spectrum and a circular orbit. Since a\ has a va lue of 



a few in the optical for the typical temperatures of He WDs (e.g., iKilic et al.ll2010al ). the 
equation above shows that the beaming amplitude is expected to be at the few 1CT 3 level, 
and therefore detectable from the ground. We note that this is not the accuracy required 
from a single measurement but the overall relative sinusoidal variability amplitude along the 
entire orbit. 

Fig.[T]shows the expected SDSS-g' beaming amplitude versus orbital period, for NLTT 1174£ 
like systems (solid line). Also plotted (dashed line) is the expected amplitude for a system 
with an inclination angle of i = 60 degrees. The figure shows that systems similar to 
NLTT 11748, but with longer periods and/or lower inclinations are also expected to show a 
beaming amplitude above the 10~ 3 level. 



3. Observations and Data Analysis 

We observed NLTT 11748 for 3.5-4.0 hr during each of the three nights of 2010 February 
7, 9 & 12 UT, obtaining data over all orbital phases. Most orbital phases were observed 
during two different nights. Our observations were done at the LCOGlEI 2.0m Faulkes 
Telescope North (FTN), located at the Haleakala Observatory, Hawaii. We used the Spectral 
Instruments 600 Series camera and the SDSS-g' filter. The camera has a 4K x 4K back- 
illuminated Fairchild Imaging CCD with a 10.' 5 x 10.'5 field of view (FOV). We used the 
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default 2x2 pixel binning mode, with an effective pixel scale of 0."304 pixel . The telescope 



pointing was adjusted to have the guiding camera FOV include a suitable guide star. Our 
observations were done while defocusing the telescope, spreading the stellar PSF over more 
pixels and increasing the overall open shutter time relative to the camera dead time, mostly 
CCD readout time. The exposure time was 120 s, resulting in a median cycle time of 143 s. 

Data taken during each night was processed the same way. Preliminary reduction, 
including bias and flat-field corrections was done using standard IRAF routines. We applied 
aperture photometry using IRAF/phot by using a few aperture sizes and eventually choosing 
a radius of 10 pixels (3."04), resulting in the smallest scatter in the final light curves of the 
target and similar brightness stars, within 1 mag from that of the target. The target's light 
curve was calibrated by dividing it by a comparison light curve — a weighted average light 
curve of four comparison stars showing no variability, whether inter- or intra-night. This 
calibration included propagation of the errors of both the target and comparison light curves. 

Out of the total of 286 exposure s seven were taken completely or partially during the 



clear outliers in the target's light curve were visually verified to result from a cosmic ray 
hitting within the photometric aperture. Therefore, a total of 277 exposures were used in 
our analysis. 

The error estimates of the calibrated light curve were directly based on the results of 
IRAF/phot. We rescaled these errors to better reflect the true noise in the data, by fitting 
a 3rd degree polynomial to the flux ratio light curves of each night, and inflating the error 
bars to make the reduced x 2 equal one. The derived scaling factors were 1.83, 1.37 & 1.53 
for the three nights, Feb. 7, 9 & 12, respectively. We note that this polynomial fit was done 
only for the purpose of rescaling the error bars. 

The target's airmass varied from «1 to ~2 during our observations, which is expected to 
affect the resulting light curve, depending on the difference in color between the target and 
comparison light curves. We removed these trends by taking the data from all three nights 
together and fitting a 2nd degree polynomial to the target's relative flux versus airmass, and 
dividing by the fit. This step resulted in the final light curve, shown in Fig. [2J 

In order to look for the beaming effect in the final light curve we simultaneously fitted 
for the three known effects (beaming, reflection and ellipsoidal). In our modeling we took 
the beaming effect to be a sine modulation at the orbital period. The reflection effect has 
the same period but is lagging in phase by tt/2, hence it is a cosine modulation at the orbital 
period. The ellipsoidal effect has half the orbital period but the same phase as the reflection 
effect, so it is a cosine modulation at the first harmonic of the period. We used the ephemeris 



3 min eclipses 




analysis. Two additional 
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of ISteinfadt et al.l ( 2010T) and fitted the phase fo lded light curve with a 5-parameter model 
( jSirko fc Paczyriskil l2003t iMazeh &: Faiglerl l2010f ) , including a sine and cosine terms for the 
orbital period and its first harmonic: 



2tt 



P 



fit) = a + aic cos I —t + ai s sin I —t + a 2c cos 



2tt 



P 



2vr 
P/2' 



t) + a 2s sin 



2vr 
P/2 



(5) 



where t = t— To, time subtracted the mid primary eclipse time. The ephemeris of lSteinfadt et al 
( 120101 ) is based on observations done during 1.5 months prior to the observations presented 
here. This resulted in a negligible uncertainty in the phase during our observations, smaller 
than 20 seconds, or of the period, justifying keeping the phase fixed in our fit. 



4. Results 

Results for the fitted coefficients are listed in the middle column of Table [TJ The fitted 
beaming relative amplitude is a\ s = (3.0 ± 0.4) x 10~ 3 , showing that the beaming signal is 
detected with a significance of > 7a. The other three fitted coefficients were found to be 
consistent with zero. The x 2 statistics of the fitted model was 283.9, and the reduced x 2 
was Xred = X 2 /d°f = 1-04 for the 272 degrees of freedom, indicating a good estimate of the 
errors. Fig. |3] presents the phase folded light curve and the fitted model is overplotted. 



4.1. Further Tests 

We have tested our results in order to verify they are not a statistical fluctuation, and 
that the beaming signal was not accidently injected into the data by our analysis. 

Assuming the light curve is non-variable (equivalent to fitting a constant value equal 
to the weighted average) results in x 2 — 351.5 for 276 degrees of freedom, and x 2 e d = 1-27, 
clearly rejecting the no variability hypothesis. A bootstrap analysis, done by randomly 
permuting the light curve points in time showed that the probability of the measured beaming 
signal being a random one is smaller than 10~ 6 . Both these tests show, in addition to the 
amplitude's high S/N, that the detected variability has a high statistical significance. 

It could be that the way we accounted for the airmass induced variation has injected 
the sine signal into the data. We tested this in two ways. First, we changed the polynomial 
degree used in fitting the relative flux versus airmass. Using a degree of 1 gave a 0.5<r increase 
in the fitted beaming amplitude and increasing the polynomial degree up to 6 decreased the 
amplitude by < 0.5a. In all these fits the other three fitted amplitudes were found to be 
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consistent with zero. Second, we examined the light curves of similar color stars in the field, 
with a B — V term within 0.1 mag from that of NLTT 11748. For those stars Xred was 
similar at the 0.01 level to that of a non- variable model, indicating they show no significant 
variability, and that the way we accounted for the airmass variation during our observations 
did not induce the detected variability. 

Another possibility is that the detected signal originated from a variability in one of the 
four comparison stars. We checked this in two ways. First, we repeated our analysis four 
times while excluding one of the four comparison stars, using only the other three. This 
changed the fitted coefficients by less than la and resulted in similar x 2 values. Second, we 
applied the same analysis for each comparison star, using only the other three as comparisons. 
As in the case of the similar color stars described above, for all four comparison stars this 
resulted in a similar \ 2 when fitting the 5-parameter model and a non-variable model, and the 
fitted beaming amplitude was consistent with zero. These two tests show that the variability 
identified in the target did not originate from the comparison stars. 

We repeated our analysis three times using data from only two of the three nights. We 
were able to detect the beaming amplitude in each of the three analyzes, although with a 
lower significance, but still higher than 4a. In each of those three analyzes the other three 
coefficients were found to be consistent with zero. 

Finally, we fitted the data also with a simpler, 2-parameter model, including a zero point 
and only the sine component corresponding to the beaming amplitude (a ls ). The result was 
the same as in the original 5-parameter model, and is listed on the second column from the 
right in Table [TJ 

4.2. Comparison with expectations 

For estimating the expected amplitude of the beaming effect in our data we first cal- 
culated it as a function of wavelength assuming a blackbody spectrum, using Eq. Q. This 
function is plotted as a solid black line in Fig. HJ We then weighted it by the SDSS-g' filter 
transmission curve and the CCD QE. The final estimates for each band are marked as filled 
circles in Fig. H] where the error bars due to the uncertainties in T e g t x and K\ are comparable 
to the size of the markers. Fig. H] also justifies our choice to observe in SDSS-g', as the 
beaming amplitude increases towards the blue, but observations in SDSS-w' are difficult at 
FTN. The expected SDSS-g' amplitude is (3.3 ± 0.1) x 10~ 3 , which is within la from the 
measured amplitude of (3.0 ± 0.4) x 10~ 3 , marked by a diamond in Fig. HI Therefore, the 
measured amplitude for the beaming effect is consistent with the expected value assuming a 
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blackbody spectrum. 



Another way to estimate the expected beaming amplitude, in a more direct manner 
and without the need to assume a blackbody spectrum is to integrate over a measured 
spectrum of the system using Eq. 2 of lyan Kerkwiik et alj (120101) . For that end we used the 
NTT/EFOSC2 spectrum obtained by iKawka fc Vennesl (120091 ). which covers the SDSS-g' 
band. This resulted in an expected amplitude of 3.7 x 10 -3 , less than 2a away from the 
measured amplitude. However, we caution that the uncertainty on this result is difficult 
to estimate. This low resolution spectrum was derived from two 20 min exposures, and 
was affected by the Doppler shift variation during the exposure, and by the RV difference 
between the two exposures, which could have been at the ~50 km s -1 level, depending on 
the exact orbital phases they were taken at. 



We estimate the amplitude of the ellipsoidal effect following iMorris &: Naftilanl (119931 ) 
to be about ^ (^) 3 , where a is the orbital semi-major axis, r-\ the radius of the primary. 



and we have ignored an order to unity coefficient (see also lLoeb fc Gaudill2003t IZucker et al. 



20071 : iMazeh fc Faiglerlboioh . This expression is < 10~ 4 for NLTT 11748. Although it is 
only a rough estimate and a more accurate one can be derived for this system, we use it only 
to show that the ellipsoidal effect is too small to be detected in our data. For the reflection 
effect we use (— ) as an upper limit, which is 4 x 10~ 5 for NLTT 11748. Therefore, both 
the ellipsoidal and reflection effects are not expected to be detected from the ground, in 
agreement with our results. 



5. Summary and Conclusions 

We have presented a ground-based measurement, carried out with a 2. 0m telescope, of 



the re lativistic beaming effect in the high mass ratio WD binary NLTT 11748 (jSteinfadt et al. 
20101 ). We showed that the measured amplitude, of (3.0 ± 0.4) x 10~ 3 , is within la from an 



approximated theoretical estimate, and within 2a from a more direct estimate using a low 
resolution spectrum of the system. The ellipsoidal and reflection effects were not detected 
in our data, in agreement with theoretical expectations. To our knowledge this is the first 
ground-based measurement of the relativistic beaming effect. 

Our measurement is a first step in testing the relation between the beaming and RV 
amplitudes in systems similar to NLTT 11748, and it suggests it can be approximated by 
assuming a blackbody spectrum. Although we were not able to obtain a high precision 
measurement of the beaming amplitude, it is identified at a high significance, of more than 
la. 
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Although short period detached He + C/O WD binar ies are rare, their discovery rate 



is increasing. Five such systems were recently discovered (jMullally et al.l 120091 ; iKilic et al. 



2010al ; iMarsh et al.l |2010| ; iKulkarni fe van Kerkwijkl 120101 ) by RV monitoring, and are ex- 
pected to show beaming amplitudes of > 2 x 10~ 3 in the optical. All five have orbital 
periods of only a few hours, shorter than that of NLTT 11748, so they can be used to fur- 
ther test the relation between the RV amplitude and the beaming amplitude, and also the 
wavelength dependency of the latter, shown in Fig. [4] for NLTT 11748. 

The evolution of a low- mass He WD requires it to have a binary companion (e.g., 
Marsh et al.lll995l ). Therefore, once a He WD is identified the beaming effect can be used 
to look for a binary companion photometrically, in a very similar way to spectroscopic 
monitoring. Since in this case the light curve becomes analogous to the RV curve, high 
precision photometry can in principle be used to model the orbit and measure the mass ratio 
of the two components. 
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Table 1. Fitted coefficients 



Coefficient 


Effect 


5-parameter model 
Fitted value 
x 10 3 


2-parameter model 
Fitted value 
x 10 3 


Expected value 
x 10 3 




Beaming 


3.0 ±0.4 


3.0 ±0.4 


3.3 ±0.1 


die 


Reffection 


0.3 ±0.4 




< 0.1 






0.1 ±0.4 






a-2c 


Ellipsoidal 


0.1 ±0.4 




< 0.1 



i=90° 

i=60° 




Period [hours] 

Fig. 1. — Expected beaming amplitude according to Eq. (jl]), for an NLTT 11748-like system 
as a function of orbital period in the SDSS-g' band. The plot was produced using the masses 
of the two WDs in NLTT 11748 and the effective temperature of the primary. The solid line 
corresponds to a completely edge-on system, close to NLTT 11748 inclination of 89.9 ± 0.1 
degrees, while the dashed line corresponds to a similar system with 2=60 degrees. The black 
filled circle marks the expected NLTT 11748 position. The plot clearly shows that similar 
systems to NLTT 11748, with a longer period and a lower inclination are also expected to 
show a beaming amplitude above the 10 -3 level, marked with a horizontal dotted line. 
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34.7 34.8 34.9 36.7 36.8 36.9 39.7 39.8 39.9 

BJD - 2455200 

Fig. 2. — Nightly light curves of NLTT 11748, observed at FTN on the nights of (from left to 
right) 2010 February 7, 9 & 12 UT. The overplotted solid line is the fitted model. Although 
the scatter in the data is relatively large, the small variability can be seen in the data. 
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Fig. 3. — NLTT 11748 light curve phase folded with the known ephemeris. Phase zero 
is time of primary eclipse. Data from different nights are marked by different colors (top 
light curve). The bottom light curve was generated by binning the points in phase, using 
0.04 phase bins, and was shifted downward for clarity The binned light curve is shown for 
visualization only, and was not used in our analysis. Overplotted solid line (black) is the 
fitted 5-parameter model, constructed from the coefficients listed in Table [TJ 
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Fig. 4. — Amplitude of the expected relativistic beaming effect for NLTT 11748 as a function 
of wavelength (black solid line), assuming a blackbody spectrum. The filled circles are the 
calculated amplitudes in each filter, accounting for the CCD QE and filter transmission 
curves. The plotted colored curves account for both the CCD QE and filter transmission. 
The SDSS-g' beaming amplitude measured here is marked by a diamond, along with error 
bars (green). 



